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Abstract . We consider the interaction of gravity, as expressed by Einstein's Equations 
of General Relativity, to other force fields. We describe some recent results, discussing 
both the mathematics, and the physical interpretations. These results concern both 
elementary particles, as well as cosmological models. (This paper describes joint work 
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1 Introduction 



We wish to describe how General Relativity modifies classical physics on two different 
scales: the scale of elementary particles, and on cosmological scales, the large scale 
structure of the Universe. We shall discuss a few different scenarios, describing the main 
ideas. The more detailed mathematical discussions can be found in our original papers. 

The mathematics of General Relativity (GR) has several features which are extremely 
interesting, and are quite different from the classical Newtonian theory of gravity. We 
will show how rich a subject GR is, both mathematically, and in its physical applications 
and interpretations. Indeed, GR is most interesting from a physical viewpoint because 
it gives a different understanding of physical phenomena. 

The main feature of Einstein's Theory of GR is that the gravitational field is the metric 
in 4-dimensional spacetime. Thus GR is actually a theory of spacetime in the sense that 
the underlying spacetime continuum is not fixed, but is allowed to vary; indeed, massive 
bodies modify the curvature of their surrounding spacetime geometry, and free particles 
traverse along geodesies determined by the metric. The second important feature of GR 
is the existence of black holes; this notion has no classical analogue. We shall see the 
role played by black holes in various aspects of this paper. Finally, we remark that GR 
in and of itself, is a very beautiful subject. This was already recognized by Einstein 
himself who in 1915, in his presentation of the theory of GR to the Prussian Academy in 
Berlin stated, "Hardly anyone who has grasped the theory, will be able to escape from 
its magic." This is undoubtedly still true today. 

2 Background Material 

In this section we shall review some basic ideas in GR, Yang-Mills (YM) equations 
and the Dirac (D) equation in four-dimensional Lorentzian spacetime. 
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A. General Relativity 



General Relativity is Einstein's Theory of Gravity, and is based on the following 3 
hypotheses. 

(Ei) The gravitational field is the metric g^ in 4-dimensional spacetime. The metric is 
assumed to be symmetric: = g^ i,j — — 3. 

(E 2 ) At each point, g^ can be diagonalized as 

gn = diag(-l, 1,1,1). 

(E 3 ) The equations should be independent of the particular coordinates involved, and 
thus the equations should be tensor equations. 

(Ei) is Einstein's brilliant insight, whereby he "geometrizes" the gravitational field, 
replacing one Newtonian potential for the gravitational field by the ten metric potential 
functions, gij, % < j. (E2) means that Special Relativity is included in GR, (and also 
that the metric tensor is everywhere non-singular), while (£3) asserts that coordinates 
are merely an artifact, and physics shouldn't depend on the choice of coordinates. 

The metric g^ = gij(x), i,j = — 3, x = (x°, x 1 , x 2 , x 3 ), x° = ct (c = speed of light), 
is a tensor defined on 4-d spacetime. Einstein's equations are ten (tensor) equations for 
the unknown metric g^ (gravitational field), and take the form 

R ij ~ \ R 9ij = vTij . (2.1) 

The left-hand side = — \ Rgij is the Einstein tensor and is a geometric quan- 
tity, depending only on g^ and its derivatives, while Tij, the energy-momentum tensor, 
represents the source of the gravitational field, and encodes the distribution of matter. 
The word "matter" in GR refers to everything which can produce a gravitational field, 
including elementary particles and electromagnetic fields. Since the Einstein tensor iden- 
tically satisfies the equation G 1 ^ = (the covariant divergence vanishes) it follows that 
on solutions of (1), Tj ;i = 0, and this in turn expresses the laws of conservation and 
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momentum, [23]. The quantities which make up the Einstein tensor Gij are given as 
follows. First from the metric tensor g^, we form the Levi-Civita connection (Christofffel 
symbols) Y\-\ 

where [g ke ] = [gke] -1 , and summation convention is employed; namely every up-down 
index is to be summed from - 3.Q Using these r&-, we construct the Riemann curvature 
tensor R} qM : 

pi _ it __ gfc , -pi pP _ pi pP 

i M ~ dx k dx e P k i £ p £ i k ' 

Finally, we can explain the terms Rij and R in G^; namely Rij = Rf s j is the Ricci tensor 
and R = g^Rij is the scalar curvature. 

The quantity a is a universal constant defined by 

87TK 

G = — 

where k is Newton's gravitational constant, and c is the speed of light (both in "appro- 
priate" units; we shall often choose units in which k — 1 = c). 

From these definitions, one immediately sees the enormous complexity of the Einstein 
equations (2.1). For this reason, we shall often seek solutions which have a high degree 
of symmetry, the aim being to make the resulting equations mathematically tractable. 

B. Black Hole Solutions 

Consider the gravitational field outside of a ball of mass M in IR 3 . Solving Einstein's 
equations Gy = 0, gives the celebrated Schwarzs child solution (1916): 



1 In a letter to a friend Einstein wrote, "I have made a great discovery in mathematics; I have 
suppressed the summation sign every time the summation must be made over an index that occurs twice 
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ds 2 = - (l - ^ c 2 dt 2 + (l - ^ 1 rfr 2 + rW, (2.2) 

where m = and <if2 2 = <i# 2 + sin 2 6><i0 2 is the standard metric on the unit 2-sphere. 
Since 2m has dimensions of length, it is called the Schwarzschild radius. Observe that 
on the sphere r = 2m, the metric becomes singular; indeed g tt = and g rr = oo. By 
transforming the metric (2.2) to so-called Kruskal coordinates, [1], one sees that the 
Schwarzchild sphere r = 2m, has the physical characteristics of a black hole: light and 
nearby particles can enter, nothing can escape, and there is an intrinsic (non-removable) 
singularity at the center r = 0. 

More generally, we could consider a metric of the form 

ds 2 = -T(r) 2 dt 2 + A-\r)dr 2 + r 2 dVl 2 , (2.3) 

where 

A( r ) = 1-^^-, 2m(r) = r(l - A(r)); 
r 

m(r) is called the "mass" function. For this metric, we define a black hole solution to be 
a solution of Einstein's equations which satisfies 

A(p) = 0, for some p > 0, 

and A{r) > if r > p. p is the radius of the black hole, and is often referred to as the 
event horizon. (Defined by the condition that it be the largest zero of the dr 2 term in 
the metric.) 

C. Yang-Mills Equations 

The Yang-Mills equations are a generalization of Maxwell's equations of electromag- 
netism. To see this, we first write Maxwell's equations in an invariant way. For this, let 
A be the scalar-valued 1-form 

A = Aidx 1 , Ai e IR , 
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called the electromagnetic potential (by physicists), or a connection (by geometers). Then, 
F, the elecromagnetic field (curvature) is the 2-form defined by 



F — dA. 



(2.4) 



In coordinates, F can be written as 



F = Fijdx 1 A dx 3 , 



dAi dAj 



dx 3 dx % 



In this set-up, Maxwell's equations take the form 



dF = , 



d*F = 0. 



The first equation follows trivially from (2.4), while for the second equation, the "star" 
is called the Hodge star operator and in four dimensions, it maps 2-forms to 2-forms and 
is defined by 



where g = det(gij), and is the completely anti-symmetric symbol defined as = 
sgn(i,j, k,£). Here, as usual, we always use the metric to raise (or lower) indices, so that 



It is important to notice that *F depends on metric gij. 

The Yang-Mills equations are a set of equations which generalize Maxwell's equations. 
Thus, to each Yang-Mills equation is associated a Lie group G (of symmetries), called 
the gauge group. For such G, we consider its Lie algebra $j, defined as being the tangent 
space at the identity of G. Now suppose that A is a 0-valued 1-form; i.e. 



where each is in q. In this (slightly!) more general case, the curvature 2-form F is 
defined by 




p l 3 




F = dA + A A A, 
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and in components, 



Notice that the commutator, [Ai, Aj] = if G is an abelian group, but is generally 
non-zero if G is a matrix group. The Yang-Mills equations are again given by 

dF = and d*F = 0. (2.5) 

These generalize Maxwell's equation because for Maxwell's equation, G = U(l), the circle 
group (U(l) = {e ld : 9 e H}), so q consists of scalars and the commutator term vanishes. 
Note that for Maxwell's equations d*F = is a linear equation for the F^s, or for the 
unknown "connection coefficients" Ai. However if G is non-abelian, say G = SU(2), 
then the Yang-Mills equations d*F = become non-linear equations for the (unknown) 
matrices Ai. 

D. The Dirac Equation in Curved Spacetime 

The Dirac equation brings in Quantum Mechanics, and is a generalization of Schrodinger 
equation to the relativistic case. It also describes the intrinsic "spin" of certain elemen- 
tary particles. The Dirac equation takes the form 

(G-m)* = 0, (2.6) 

where G is the Dirac operator and ^ is a complex- valued 4- vector called wave the function 
(spinor) of a fermionF] (Dirac particle) of mass m. The Dirac operator G can be written 

as 

G = i&(x)-^+B(x) 

where the G J are 4x4 matrices called the Dirac matrices and B is a 4 x 4 matrix. The 
Dirac matrices G^ and the Lorentzian metric (defined on 4-d spacetime) are related 
by 



gi"l=±{G>,G k }, (2.7) 



2 Fermions are distinguished from bosons in that fermions have half-integral multiples of spin, while 
bosons have integral multiples of spin. Only fermions obey the Pauli Exclusion Principle. Examples of 
fermions are protons and antiprotons, neutrinos, electrons and positrons, quarks, and leptons. Examples 
of bosons are mesons, photons, and pions. 
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where \ G j , G k \ is the anti- commutator 



\G 3 ,G k } = G j G k + G k G j . 

The Dirac matrices thus depend on the metric. 

Now let H be any space-like hypersurface in IR 4 , with future directed normal vector 
field v = (z/j), and let dfi be the invariant measure on H induced by the metric gij. Define 
a scalar product on solutions, $ of the Dirac equation by 

= f ^/G^Ujdfi . 



This scalar product is positive definite, and, as a consequence of current conservation 
(c.f. [4]) 

V^G j $ = 0, 

it is independent of the choice of the hypersurface H. By direct generalization of the 
expression 

\j>y\j/ — |\i/| 2 

in flat Minkowski space given originally by Dirac, ([2]), where 

1 
-1 

in our case \1> is the adjoint spinor, defined by 

1 
-1 

1 is the 2x2 identity matrix, and ^* denotes complex conjugation, 

*G%i/j (2.8) 
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is interpreted to be the probability density of the Dirac particle. We normalize solutions 
of the Dirac equation by requiring 

= 1. (2.9) 
If the Lorentzian metric admits a black-hole solution at r = p > 0, (2.9) is replaced 

by 

^G^Vjdfi < oo (2.10) 

{t=const., r>ro} 

for all r > p. For normalized solutions of the Dirac equation, the integral (2.10) gives 
the probability for the Dirac particle to be in the region r > r . We shall only consider 
Dirac particles which lie outside the event horizon. This is because the probability 
density (2.8) is not necessarily positive inside the event horizon so that a positive infinite 
contribution of the integral outside the event horizon can be compensated by a negative 
infinite contribution inside the event horizon. We thus demand that the integral (2.10) 
away from, and outside of the event horizon, be finite. 

In the paper [4] , it is shown that the Dirac matrices G j can be chosen to be any 4x4 
matrices which are Hermitian with respect to the scalar product 

\]>|cl>\ = / (I^a/I^I dx 
1 Jn 4 

and satisfy the anti-commutation relations (2.7). This gives us more flexibility in choosing 
the Dirac matrices. Namely, the relations (2.7) do not uniquely determine the Dirac 
matrices in curved spacetime. However it was proved in [4] that all choices of Dirac 
matrices satisfying (2.7) yield unitarily equivalent Dirac operators. 

With this introduction to the Einstein equations, the YM equations and Dirac's equa- 
tion, we turn to our first topic, the coupled Einstein-Dirac- Yang/Mills (EDYM) equa- 
tions. 
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3 Einstein-Dirac- Yang/Mills Equations 

(joint work with F. Finster and S.T.-Yau) 



The EDYM equations are obtained by varying the action 



S 



167TK 



R + V(G -m)tf 



1 



167re" 



g\ dx 



over Lorentzian metrics g^, g- valued Yang/Mills connections Aidx 1 (where g is the Lie 
algebra of the given gauge group), and 4-spinors x &. Here e is the so-called Yang/Mills 
coupling constant, and k is Newton's gravitational constant. In general, the EDYM equa- 
tions are extremely complicated. We specialize here to spherically symmetric solutions 
depending only on the radius r = \x\. Taking the gauge group as SU (2)0 together with a 
special ansatz for the spinors, we can simplify the Dirac equations from complex 4-spinors 
to the case of an equation for real 2-spinors (a, /?)*, [8]. In this set-up the (SU(2)) EDYM 
equations take the following form: 

7 W/r -{m + uo)T\ 





1 


\f>i 


1 " 



— m + ujT 



-W/r 



a 



(3 



(3.1; 



2rA 



rA' 



T 



I — A — 



K (1 



- 2KwT 2 (a 2 + f3 2 ) - ^A(W'f (3.2) 



I + .4 + 4^ — + 2mT(a 2 - (3 2 ) - 2kujT 2 {oi 2 + (3 2 ) (3.3) 



g2 /y 2 



+ 4k— Waf3 - ^rA(W') 2 



r 2 AW" = -(1 - W 2 )W + e 2 rTa(3 



A'T - 2AV 
2T 



W. 



(3.4) 



Here m is the rest mass of the Dirac particle, u is its energy, and W is the unknown con- 
nection coefficient. Equation (3.1) is the Dirac equation, (3.2) and (3.3) are the Einstein 

3 The gauge group SU(2) corresponds to the weak nuclear force, while SU(3) corresponds to the 
strong nuclear force. 
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equations (for the Lorentzian metric of the form (2.3)), while (3.4) is the Yang/Mills 
equation. Notice that equations (3.1), (3.3) and (3.5) all become singular when A = 0; 
i.e., at a black-hole. (This is special to Einstein's equations, and is the first significant 
difference that we see from classical equations.) 

Smooth solutions of (3.1) - (3.4) are those which are defined for all r > 0, and are 
called particle-like solutions. Black hole (BH) solutions correspond to solutions having 
A(p) = for some p > 0, and are defined for all r > p, with A(r) > 0, if r > p. In 
addition solutions of (3.1)- (3.4) are required to satisfy certain initial condictions, (see [8, 
9]), together with the following global conditions : 

/ (a 2 + l3 2 )—j= dr = 1 (for particle — like solutions), (3.5) 
Jo v A 



< / (a 2 + [3 2 )—= dr < oo for each r > p, 

Jro VA 

(for BH solutions), 



(3.6) 



lim 



-DC' 



lim r(l — A(r)) < oo , (3.7) 

lim T(r) = 1 (3.8) 

r— >oo 

(w(r),W"(r))), is finite. (3.9) 



Conditions (3.5) and (3.6) state that the spinors (wave functions) are normalizable, (3.7) 
means that the fermions have finite (ADM) mass, (3.8) (together with (3.7)) means that 
the gravitational field is asymptotically Minkowskian, while (3.9) asserts that the YM 
field is well-behaved. We now state our main results, and discuss them afterward. 



Theorem 3.1 ([8]): There exist stable particle-like solutions of the EDYM equations for 
arbitrarily weak gravitational coupling constant: m 2 k/e 2 — > 0. 

Theorem 3.2 ([9]): Every black hole solution of the EDYM equations cannot be normal- 
ized; namely, the spinors must vanish identically outside of the black hole. 



11 



Both of these results are unexpected and rather surprising. 

Theorem 3.1 is obtained using numerical methods. This (stability) result has been 
shown to be false for the (SU(2)), Einstein- Yang/Mills (EYM) equations, [21, 22], and 
shows that Dirac particles in a gravitational field form bound states if an additional 
strong coupling to a non-abelian Yang-Mills field is taken into account. In [8], we employ 
some new numerical techniques, including a multi-parameter "shooting method" , and a 
scaling technique introduced in [5]. Our result shows that weak as gravity is (e.g., it 
is 10 21 times weaker than electromagnetism), it still has a regularizing effect. Although 
gravity is not renormalizable (which means that the problem cannot be treated in a 
perturbation expansion), our solutions of the EDYM equations are regular and well- 
behaved for arbitrarily weak gravitational coupling. Our stability result is obtained by 
using Conley Index techniques, [17]. 

The proof of Theorem 3.2 is given in [9]. The result shows that the EDYM equations 
do not exhibit normalizable BH solutions. Thus, in the presence of quantum mechanical 
Dirac particles, static, spherically symmetric BH's do not exist. Another interpretation 
of our result is that it indicates that Dirac particles can only either disappear into the 
BH or escape to infinity. We shall discuss this further in the next section. 

4 Decay Rates and Probability Estimates for 
Massive Dirac Particles in a Charged, 
Rotating Black Hole Geometry 

(joint work with F. Finster, N. Kamran, and S.-T. Yau) 

Consider the physical situation of a Dirac particle (fermion, electron, proton, etc.) of 
mass m > outside of a charged, rotating black hole. Question: What is its long time 
behavior, as t — > oo? 

In the classical case (no Dirac equation, no quantum mechanical considerations), the 
particle can remain in a time-periodic orbit around the black hole (see [3]). In this 
section we shall show that our previous "indication" , (mentioned at the end of the last 
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section) is valid. That is, if quantum mechanical effects are taken into account, the 
picture completely changes, and the Dirac particle either enters the black or hole, or 
tends to infinity; no other possibilities can occur. 



In the paper [7], we proved that there are no normalizable, time periodic solutions of 
the Dirac equation in a Reissner-Nordstrom^ black-hole background; in particular there 
are no static solutions of the Dirac equation in such a background metric. This result 
was extended in [10] to the case in which the background geometry is that of a charged, 
rotating black hole. In this section we shall discuss such black holes and we shall show 
that our above "indications" are correct; i.e., in such a background geometry, a Dirac 
particle either tends to infinity, or enters the black hole - there are no other possibilities. 
In addition we shall obtain probability estimates on which of these two possibilities 
occurs. We begin with a brief discussion of the geometry of a charged rotating black 
hole. 



In 1963, R. Kerr found a solution of Einstein's equations corresponding to a rotating 
black hole [15]. This result was generalized by Kerr and Newman; (see [11]), to a charged 
rotating black hole (called the Newman- Penrose solution of Einstein's equations): 



ds 2 
where 



A 

U 



(dt — a sin 2 



u ( ^ + de 2 



sin 



U 



adt — (r 2 + a 2 )d(f) 



(4.1) 



U{r,9) 
A(r) 



2 i 2 2/1 

r + a cos 6 
r 2 - 2Mr + a 



Q 2 



and M, aM, and Q denote the mass, angular momentum, and charge of the black hole, 
respectively; a is the angular velocity of the black hole. From our earlier remarks, the 
largest root of A corresponds to the event horizon of the black hole. We assume^ 

M 2 >a 2 + Q 2 , (4.2) 



then 



P 



a- 



Q 2 



4 The Reissner-Nordstrom solution is an extension of the Schwazschild solution, and is a static solution 
of the coupled Einstein-Maxwell equations; see [1]. 
5 This is the so-called "non-extreme" black hole. 
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defines the event horizon, the boundary of the charged, rotating, black hole (CRBH). 



Now consider the Cauchy Problem for a Dirac particle of mass m > 0, and charge q, 
in this CRBH background geometry, (see [1 1] , where we write the Dirac equation in this 
different way, taking the spin connection into account) having initial data outside of the 
event horizon r = p: 

(iy'Dj -m)V(t,x) = (4.3) 
V(0,x) = V (x), \x\> p. (4.4) 

Here is our first result. 

Theorem 4.1 : Let 5 > be given, and let R > p + 5. Consider the (annular) region in 
1R 3 

Ks, R = {p + 5<r<R}. 

Then the probability for the Dirac particle to be inside Ks,r tends to zero as t — > oo; 
namely, 

lim / (#y#)(i,xW^ = 0. (4.5) 



Thus with probability one, the Dirac particle leaves every bounded set in IR 3 . The 
proof of this theorem is based on a result of Chandrasekhar, [3] , who showed that the 
Dirac equation in the CRBH geometry can be separated into ODE's. Moreover, using 
this, we can write the Dirac propagator as (see [11, 12]) 

*(t, X ) = -J2 e ~ wt E tT^rt*) <*n*o>^, (4.6) 

71 fc,neZ^~°° a,b=l 

where the (positive) scalar product is defined in [11]. Here k and n are generalized 
angular momentum numbers (arising from Chandrasekhar's separation theorem), ^^° n 
are solutions of the Dirac equation, and the t^ n are generalized transmission coefficients. 
These functions can be written in terms of the fundamental solutions of the resulting 
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ODE's. For |u;| > m, near the event horizon, the ty^ 71 go over to spherical waves. 
corresponds to incoming waves (moving towards the black hole), and \fr% wn corresponds 
to outgoing waves (moving away from the black hole). Near infinity, the ty^n again go 
over to spherical waves. For \u\ < m, the fundamental solutions (for a = 1,2) near the 
event horizon are linear combinations of both incoming and outgoing waves. 
^2 lJjn ) decays (resp. grows) exponentially at infinity. The theorem is proved using the 
representation (4.6). 

Theorem 4.1 implies that the Dirac particle must leave every bounded set in IR 3 . 
Equivalently the Dirac wave function \1/ decays to zero in L^ c , outside and away from the 
event horizon. It follows that the Dirac particle must eventually either disappear into the 
black hole, or escape to infinity; these are the only possibilities. This raises the following 
questions: 

(A) What is the likelihood of each of these two possibilities?, 
and 

(B) What is the rate of decay of the Dirac wave function on a compact subset of IR 3 
outside of the black hole? 

We consider question (B) first. To this end, let q denote the charge of the Dirac 
particle and assume that the Dirac particle has "small charge" , in the sense thatg 

mM > \qQ\. (4.7) 

We also assume that the angular momentum of the initial data is bounded. This means 
that the summation in (4.6) is over a bounded set of integers n and k] say 

\k\ < ko and |n| < uq. 

6 This means that the gravitational attraction is the dominant force far from the black hole. Indeed, 
at large distances from the black hole the metric is asymptotically Minkowskian, and the gravitational 
and electromagnetic fields can be described by the Newtonian limit. In this limit the gravitational and 
electrical forces are mM/r 2 and qQ/r 2 , the Newton and Coulomb laws, respectively. 



15 



We then have the following theorem; (see [12] for the proof). 



Theorem 4.2 : (Decay Rates): Consider the Cauchy problem (4-3), (4-4) f or the Dirac 
equation in the non-extreme CRBH background geometry, with small charge. Assume 
that the Cauchy data is smooth with compact support outside of the event horizon, 
and has bounded angular momentum. 

(i) If for any k and n 

limsup |(^r n |*o}| ^ 0, or liminf \(^ n | tf )| ^ 0, (4.8) 

then 

\V(t,x)\ =cr 5/6 + 0(r 5/6 ~ £ ) as t^oo, (4.9) 
where c = c(x) ^ and any e < ^. 

(ii) If for all k, n, and a = 1, 2, (^^ n | \$/ Q ^ — q f or a \\ u in a neighborhood of ± m, 
then for any fixed x, ^/(x,t) decays rapidly in t, (faster than any polynomial). 

For a quantum mechanical particle, the Heisenberg Uncertainty Principle implies 
that its kinetic energy is not precisely known. Since the associated Hamiltonian has 
continuous spectrum (the operator is defined on an unbounded interval), the energy uj 
of the Dirac particle is a continuous parameter. Condition (4.8) means that the initial 
energy distribution has outgoing components near w = mor near u = —m. Our theorem 
implies that in this case, the decay rate is t~ 5 / 6 . That is, the decay rate t~ h ^ quantifies 
the effect of the black hole's attraction, on the long time behavior of the Dirac particle. 

In flat Minkowski space solutions of the wave equation having compactly supported 
data, decay rapidly (Hughyens principle). In the CRBH geometry, the Dirac wave func- 
tion \1> behaves near infinity like solutions of the wave equation in flat Minkowski space, 
and near the black hole, it should behave like a massless particle, which decays like t~ 3 / 2 
(see [12]). Thus one would expect (by "interpolation") that the Dirac particle in the 
CRBH geometry should decay at least as fast as t~ 3 / 2 . However, Theorem 4.2 shows 
that the naive picture is incorrect, and that the gravitational field affects the behavior of 
the Dirac particle in a far more subtle way. 
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Finally, one can ask why uj = ± m play such special roles. To answer this, note that 
for a particle of mass m in flat space, its energy uj — m, (or — m since the Dirac equation 
has negative energy solutions). If the particle's momentum is non-zero, then oj > m. Thus 
in flat space, its energy is outside of the interval (— m, m). But a quantum mechanical 
particle can have a continuous energy distribution, but only outside of (— m, m) (this is 
called the "energy gap" or "mass gap"). However, in the presence of a gravitational field, 
the energy distribution can be supported on the entire real line, but the dominant term 
is near ± m. 

We now turn to question (A), and discuss the probability for our Dirac particle to 
escape to infinity, or to enter the black hole. Let R> p, (the event horizon), and define 
the probability p of the Dirac particle to escape to infinity by 

p=lim / (^y i ^)(t,x)i/ j dti. (4.10) 

Jr>R 

Note that p is independent of R; namely, if R2 > R\ > p, and 



Pi = lim / (^'j j ty)(t,x)v j dx, i = l,2, 



then 



Pi - Pi — lim / - / = lim / = 0, 



t^oo V l r>Ri J r>R J t^oo jRi 



R2 



by Theorem 4.1; thus p\ = p 2 and this proves our assertion. In [12] we prove that 

1 . _\ / . v 2 



P =- y [ ( - 2 

vr .frri ./fl\[_ m>m ] V 2 



\k\<k 
\n\<n 



■■12 



du. 



Accordingly, \ —p gives the probability that the Dirac particle enters the black hole. The 
following theorem gives conditions under which p — 0, p = 1, or < p < 1, in terms of 
the initial energy distribution. 



Theorem 4.3 ([12]): Consider the Cauchy problem as in Theorem \.2, with initial data 
normalized by ^ |^ ^ = 1. 
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(i) If the outgoing initial energy distribution satisfies ^^"l^oy 7^ for some uj with 
\cu\ > m, then p > 0. 



(ii) 


If the initial energy distribution satisfies for a = 1 or 2, ^^"l^oy 


7^ for some 


UJ, 


\uj\ < m, then p < 1. 




(iii) 


// the initial energy distribution is supported in the interval [—m, m] 


, then p = 0. 


(iv) 


If (4.8) holds, then < p < 1. 




(v) 


p = 1 if and only if for all k,u and n, the following conditions hold: 





n |^ ^ = 0, if M < m 
^f n |^ ^ = -2^ n ^^"|^ ^ if uj>m. 



Note that in case (i), the Dirac particle has a non-zero probability of escaping to 
infinity; in case (ii) the particle can enter the black hole, while in case (iii) the components 
of the initial energy distribution, which lead to the decay rate t~ 5//6 , do not have enough 
energy to allow the Dirac particle to tend to infinity. Indeed, ^r| wn for \uj\ > m is outgoing 
near the event horizon, and so the Dirac particle resists the gravitational attraction for a 
while, but it eventually gets turned around and is driven into the black hole. In case (iv), 
the particle has a positive probability of either escaping to infinity, or entering the black 
hole. Finally we can understand why p = 1 for special choices of initial data. To obtain 
such data, consider the special situation where a Dirac particle at time t = —00 comes 
in from spatial infinity. Taking \l/(0,x) as initial data and reversing the time direction, 
the solution of this Cauchy problem clearly escapes to infinity with probability 1. 
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5 Shock- Waves, Cosmology, and Black Holes 



(joint work with B. Temple) 

In the previous sections we have discussed the behavior of elementary particles in 
a gravitational field; that is, how General Relativity affects matter on small scales. We 
now change directions and discuss GR on large scales, that is we shall discuss some recent 
astrophysical implications of GR. 

We again consider the Einstein equations (2.1) where now the energy- momentum 
tensor is that of a perfect fluid; i.e., a fluid in which dissipation effects are neglected: 

T ij = ip + p)uiUj + P9ij- (5-1) 

Here p and p are the density and pressure respectively, of the fluid, and u = (u°, u 1 , u 2 , u 3 ) 
is its 4-velocity, where as usual, Ui = g^uK The Einstein equations (5.1) describe the 
simultaneous evolution of the fluid and the gravitational field: matter is the source of 
spacetime curvature in Einstein's theory. 

Shock-waves are relevant here since the Relativistic Euler Equations are "embedded" 
within the Einstein equations Gij = aTy. This is because the Bianchi Identities in 
geometry, (cf. [23]), imply that the (covariant) divergence of the Einstein tensor G^, 
vanishes identically: Div = 0, so Div = 0, and this latter equation is precisely the 
Relativistic Euler Equations, expressing the conservation of energy and momentum. 

We pause to make the following remarks: 

(A) There is no "Glimm's Theorem" ([13, 17]) in GR, not yet anyway. However an 
important first step in this direction was recently made by Groah and Temple, see 
[14]. 

(B) No proof is known which demonstrates that in GR, shock-waves form from smooth 
data (see [16, 17]). 

(C) In GR, the initial data cannot be arbitrarily prescribed on a non-characteristic 
surface, and must satisfy additional constraints which are imposed by geometrical 
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considerations (Bianchi Identities). These constraints take the form of a coupled 
set of 4 nonlinear elliptic equations. There are deep unresolved issues concerning 
these constraint equations. 



In [18, 19], we constructed the first examples of shock-wave solutions in GR. The 
shock- waves are in the fluid variables and in the metric (gravitational field). This was 
done by matching two different spherically symmetric metrics Lipschitz continously across 
a (spherically symmetric) surface of discontinuity for the fluid variables (shock- wave). 
The Inner Metric is the Friedmann-Robertson- Walker (FRW) metric of cosmology. Its 
line element is given by 

dr 2 



ds 2 = -dt 2 + R(ty 



1 — kr 2 



r 2 (d# 2 + sin 2 6d(f) 2 



(5.2) 



where by a suitable re-scaling of r, we can take k = 0, 1, or — 1. This metric is clearly 
spherically symmetric, and describes a homogeneous, isotropic spacetime (no preferred 
point or direction). It models an expanding universe. There is a singularity (R = 0) 
in backwards time (i.e., earlier than present time) which corresponds to the Big-Bang. 
The function R(t) determines the red-shift factor for distant astronomical objects. (The 
red-shift is what astronomers see and can measure. It is actually the Doppler effect 
applied to astrophysics). The red shift gives information on distance, mass, and chemical 
composition of the objects; see [23] for a discussion of these things. 

The Outer Metric is the Tolman-Oppenheimer-Volkoff (TOV) metric, 

d s 2 = -B(r)dP ^P^j 1 df 2 + f 2 (d6 2 + sin 2 Od^ . (5.3) 

It is both time-independent and spherically symmetric, and models a static spacetime 
(or the interior of a star). The function M(f) denotes the mass inside a ball of radius f, 
and is given by 



rr 

M(f) = / Ans 2 p(s)ds, 
Jo 



where p is the density function. Both the FRW and TOV metrics satisfy Einstein's 
equations for a perfect fluid; see [23]. 



20 



For cosmology, we assume that we are given an FRW metric (which models well our 
expanding Universe), and we solve the differential equations for an unknown TOV metric 
which matches the given FRW metric across an expanding (outgoing) shock wave for the 
fluid variables. Our idea is to match the given {k = 0) FRW metric with equation of 
state given by 

Pmatter 0, Pradiation oX" , 

where a is the Stefan-Boltzmann constant, and T denotes the temperature.^ In our 
shock-wave model for cosmology, this model is not an expansion of the entire universe, 
but rather a limited expansion into an outer ambient static spacetime modeled by the 
TOV metric. 

In [20] , we constructed an explicit solution of this problem, and we showed that there 
are unique density and pressure profiles for each initial radiation density. In our model, 
the initial Big-Bang is more like a classical fluid-dynamical shock-wave explosion. The 
mathematical theory of shock- waves, [17] implies that many solutions decay as t — > oo 
to the same shock-wave. Thus in our shock-wave cosmology, we cannot recover all the 
information about the early details of the explosion, (due to entropy increase across 
shock- waves). 

However, there is a problem with our shock-wave cosmology model; namely our shock- 
wave is not sufficiently far enough out. Indeed, in [20], we calculated the present shock 
position to be approximately at distance (Hubble length) 

— « 10 10 light years 
Ho 

away, where H(t) = R/R(t) is the Hubble "constant", and H = H (present time). But 
we observe quasars and distant galaxies at distances of 10 11 — 10 12 light years away. For 
a long time we were puzzled as to why we cannot get our shock- wave further out. But 
we have recently answered this question in a most unexpected way. Indeed, we recently 
proved that if the shock wave is beyond the critical length scale 

jjjj^ (Hubble length) 
7 This equation of state is generally taken to be the equation of state for the Universe at present time. 
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for the FRW metric at any time t > 0, then the Universe lies "inside a black hole", in the 
sense that (1 — ^) changes sign at the critical length scale for the FRW metric. That is, 
shock matching outside a black hole can only succeed for shocks lying inside one Hubble 
length from the center of the FRW metric at any fixed time. This is a rather surprising 
connection between shock-wave cosmology and black holes. 

We also have proved that radiation admitted from a galaxy (or star) at the instant 
when it lies a distance of exactly one Hubble length from the FRW origin will be observed 
at the origin (at a later time) to be infinitely red shifted. Thus an observer at the FRW 
origin will see the shock wave "fading out of view" as t — > oo. 

We are in the process of modifying our shock-wave cosmology so as to allow shock 
waves beyond Hubble length. 
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